Isoscaling in Peripheral Nuclear Collisions around the Fermi Energy and a Signal of 
Chemical Separation from its Excitation Energy Dependence 
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The isoscaling is investigated using the fragment yield data from fully reconstructed quasi- 
projectiles observed in peripheral collisions of 28 Si with 124 ' 112 Sn at projectile energies 30 and 50 
MeV/nucleon. The excitation energy dependence of the isoscaling parameter is observed which 
is independent of beam energy. For a given quasi-projectile produced in reactions with different 
targets no isoscaling is observed. The isoscaling thus reflects the level of N/Z-equilibration in reac- 
tions with different targets represented by the initial quasi-projectile samples. The excitation energy 
dependence of the isoscaling parameter f3' , corrected for the trivial 1/T temperature dependence, 
does not follow the trend of the homogeneous system above 4 MeV/nucleon thus possibly signaling 
the onset of separation into isospin asymmetric dilute and isospin symmetric dense phase. 

PACS numbers: 25.70.Mn,25.70.Lm,25.70.Pq 



The isotopic composition of nuclear reaction products 
provides important information on the reaction dynam- 
ics and the possible occurrence of a phase transition in 
the asymmetric nuclear matter 0, H , which is supposed 
to lead to separation into a symmetric dense phase and 
asymmetric dilute phase. It has been discussed in lit- 
erature [J. [3. 0. El to what extent such a phase transi- 
tion is generated by fluctuations of density or concentra- 
tion, typically suggesting a coupling of different instabil- 
ity modes. The N/Z ( neutron to proton ratio ) degree 
of freedom and its equilibration was studied experimen- 
tally in detailed measurements of the isotopic distribu- 
tions of emitted fragments 0, H U El El El El Iso- 
topically resolved data in the region Z—2-8 have revealed 
systematic trends, which however are substantially af- 
fected by the decay of the excited primary fragments. 
It has recently been shown EH that the effect of se- 
quential decay of primary fragments can be overcome 
by comparing the yields of fragments from two similar 
reactions. For statistical fragment production in two re- 
actions with different isospin asymmetry, but at the same 
temperature, the ratio R,2i(N,Z) of the yields of a given 
fragment (N,Z) exhibits an exponential dependence on 
N and Z. This scaling behavior is termed isoscaling E3 
and has been observed in a variety of reactions under the 
conditions of statistical emission and equal temperature 

El El El El 

In this study, we present an isosca ling analysis of a 
recent fragment yield data 0, 0, Il9l |2C| obtained 
by charged-particle calorimetry of hot thermally equili- 
brated quasi-projectiles from the reactions 28 Si+ 124,112 Sn 
at 30 and 50 MeV/nucleon. The angular coverage 
and granularity of the forward-angle multidetector ar- 
ray FAUST |2J allowed full reconstruction of the hot 
quasiprojectiles from mid-peripheral binary collisions 
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El- Excellent description of the fragment observables 
was obtained using the model of deep-inelastic transfer 
( DIT ) [2I for the early stage of collisions and the sta- 
tistical multifragmentation model ( SMM ) 23] for de- 
excitation. The large N/Z range of quasi-projectiles al- 
lowed unique studies of the dependence of fragment ob- 
servables on quasi-projectile N/Z. An increasingly inho- 
mogeneous N/Z-distribution between light charged parti- 
cles and heavier fragments was observed with increasing 
isospin asymmetry [19| ■ An exponential scaling of the iso- 
baric ratio Y( 3 H)/Y( 3 He) with the quasi-projectile N/Z 
was observed and the temperature was extracted from 
the logarithmic slope [20|. In the present work we ex- 
amine the isoscaling phenomenon in unique conditions 
of highly-selective data. Typically, the investigations of 
isoscaling focused on yields of light fragments with Z=2- 
8 originating from de-excitation of massive hot systems 
produced using reactions of mass sym metric projectile 
and target at intermediate energies [lilialig or by re- 
actions of high-energy light particle with massive target 
nucleus El • I n a recent article El j we report the inves- 
tigations of isoscaling using the heavy residue data from 
the reactions of 25 MeV/nucleon 86 Kr projectiles with 
124 Sn, 112 Sn and 64 Ni, 58 Ni targets. Here we present the 
investigation of the isoscaling phenomenon on the full 
sample of fragments emitted by the hot thermally equi- 
librated quasi-projectiles with mass A=20-30. 

A thermally equilibrated system undergoing statistical 
decay can be, within grand-canonical approach, charac- 
terized by a yield of fragments with neutron and proton 
numbers N and Z I2II25II: 



Y{N, Z) = F(N, Z) exp[B(N, Z)/T] exp(A^ n /T+Z Mp /T) 

(1) 

where F(N, Z) represents contribution due to the sec- 
ondary decay from particle stable and unstable states to 
the ground state; fj, n and /x p are the free neutron and 
proton chemical potentials; B(N, Z) is the ground state 
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binding energy of the fragment, and T is the tempera- 
ture. 

The ratio of the isotope yields from two different 
systems, having similar excitation energies and similar 
masses, but differing only in N/Z, cancels out the effect of 
secondary decay and provides information about the ex- 
cited primary fragments [14j . Within the grand-canonical 
approximation ( Eq. □), the ratio Y 2 {N, Z)/Y 1 (N, Z) as- 
sumes the form 



R 21 (N, Z) = Y 2 (N, Z)/Y 1 {N, Z) = C7exp(aV + (3Z) 

(2) 

with a — A/i„/T and [3 — A/i p /T, with A[i n and Afi p 
being the differences in the free neutron and proton chem- 
ical potentials of the fragmenting systems. C is an overall 
normalization constant. Alternatively, as introduced by 
Botvina et al. 01 > the isoscaling can be expressed as 



R 21 (N, Z)=Y 2 {N, Z)/Y 1 (N,Z)=Cexp(a'A+f3'(N-Z)) 

(3) 

thus introducing the scaling parameters which can be 
related to the isoscalar and isovector components of free 
nucleon chemical potential since a' = A(fi n + /i p )/2T 
and p' — A(/i„ — fi p )/2T. In this variant of isoscaling, 
the mass dependence is rather weak and the information 
is concentrated in the isoscaling parameter /3'. 
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FIG. 1: The isoscaling plots from the reactions of 
28g i+ i24,ii2 Sn at 50 MeV/nucleon for the full set of iso- 
topically resolved quasi-projectiles and for five bins of quasi- 
projectile excitation energy per nucleon. 

In Figs, n and [21 we present the isoscaling data 
from statistical decay of hot quasiprojectiles from the 
reactions 28 Si+ 124,112 Sn at projectile energy 50 and 30 
MeV/nucleon, respectively. Observed charged particles 
with Z<5 have been isotopically resolved and the total 



observed charge was close to the charge of the projectile 
( Z=I2-15 ) 0,^3- The observed fragment data pro- 
vide full information ( with the exception of emitted neu- 
trons ) on the de-excitation of thermally equilibrated hot 
quasi-projectiles with known mass ( A=20-30 ), charge, 
velocity and excitation energy [lH . In Figs. Hland|21 the 
isoscaling plots are presented not only for the full data 
but also for five bins of excitation energy. The isoscaling 
occurs at both projectile energies, both for the inclusive 
data and for the sub-sets with different excitation ener- 
gies. The slope depends on the excitation energy in a 
similar way as the slope of the dependence of isobaric ra- 
tio Y( 3 H)/Y( 3 He) on the quasi-projectile N/Z observed 
in 20] . On the other hand, the isoscaling parameters for 
individual excitation energy bins do not depend on the 
projectile energy and the difference of the slopes of the 
inclusive data is given by the different excitation energy 
distributions of hot quasi-projectiles at two projectile en- 
ergies. 
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FIG. 2: As in Fig. [T]but for 30 MeV/nucleon. 

In Fig. |3| are presented isoscaling plots for three se- 
lected quasiprojectiles 24 ' 26 - 28 Al ( representing the cen- 
tral part of quasiprojectile mass and charge distribu- 



tion ) in the reactions of 8 Si- 



-Sn targets at pro- 



jectile energy 50 MeV/nucleon. Isoscaling behavior is 
not observed which means that no memory of the en- 
trance channel interaction leading to production of a 
given quasi-projectile remains and the isoscaling anal- 
ysis is not sensitive to the statistical fluctuations in de- 
excitation. The slight hint of isoscaling in the case of 
28 Al may be due to low statistics or due to larger uncer- 
tainty in the number of unobserved neutrons of the re- 
constructed neutron-rich quasiprojectile 28 Al caused by 
more intense neutron emission than for proton-rich iso- 
topes 24 ' 26 A1. A behavior similar in appearance to Fig. |3| 
was reported by Tsang et al. 0] using the yields of frag- 
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ments with Z=3-8 from the reaction of 16 with 197 Au, 
208 Pb targets at projectile energy 20 MeV/nucleon [2^| . 
measured at 40°. There, the absence of isoscaling be- 
havior suggests that the sample of original prefragments 
is independent of the target isospin asymmetry, what for 
such inclusive data suggests a geometric reaction scenario 
such as breakup or abrasion and the observed fragments 
are most probably the spectator parts of the original pro- 
jectile, ejected to a large laboratory angle by the recoil 
in the projectile breakup. 
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FIG. 3: The isoscaling plots from the reactions of 
28 Si+ i24,ii2g n at 5Q MeV/nucleon for three reconstructed 
quasi-projectiles. Symbols as in Figs. 11121 

As suggested by Botvina et al. \vj\. the isoscaling pa- 
rameter trends are similar to the observed double iso- 
tope ratio temperature. The temperature for the current 
data was studied in detail in [2(j and we use the results 
of that work for comparison with the isoscaling param- 
eters. In Fig. 21 we present the observable (3'T, can- 
celing out the trivial 1/T-dependence of the isoscaling 
parameter. In Fig. 2t we use the double isotope ratio 
temperature [24| from the isotopic ratios Y( 2 H)/Y( 3 H) 
and Y( 3 He)/Y( 4 He) and in Fig. 0Jj we use the tem- 
perature obtained from the slope of the dependence of 
Y( 3 H)/Y( 3 He) on the quasi-projectile N/Z ( mirror nu- 
cleus temperature ) [2(|. Grossly, the (3'T should reflect 
the level of N/Z equilibration which is known for 
the current data where full N/Z-equilibration was 
not reached. The horizontal lines represent the estimate 
of the zero temperature values of (3'T of reconstructed 
quasi-projectiles using the known N/Z equilibration [l"T| . 
We estimate the zero temperature proton and neutron 
chemical potentials by the proton and neutron separa- 
tion energies ( /Lt Pi „ ~ —S Pl „ ). Similar approximation 
was used to extract the mirror nucleus temperature |2Cj . 
The zero temperature estimates of (3'T for all excitation 
energy bins are consistent with horizontal lines which is 
a result of the selection of quasiprojectiles with Z=12-15. 
Any experimental deviations from the zero temperature 



values of (3'T can be understood as a non-trivial depen- 
dence representing the details of de-excitation. The thick 
and thin lines represent the (3'T estimates with and with- 
out the correction for neutron emission from backtracing 
of DIT/SMM simulations p). The dashed/full lines rep- 
resent the projectile energies 30 and 50 MeV/nucleon, 
respectively. The open circles/squares represent the ex- 
perimental data at 30 and 50 MeV/nucleon while the 
full squares represent the combined data from both ener- 
gies. The experimental dependences of (3'T on excitation 
energies indicate an initial decrease at low excitation en- 
ergies ( consistent with expansion of the homogeneous 
excited source ), further they exhibit a turning-point at 
4 MeV/nucleon followed by increase toward the zero tem- 
perature value at 6-7 MeV/nucleon. The trend is similar 
at both projectile energies independent of the method 
of temperature extraction and quality of the statistical 
sample. The existence of such a turning-point can be 
understood in the context of work |20|. where the via- 
bility of the approximation /i Pi „ ps —S Pjn at high tem- 
peratures was supported by a possible counterbalance of 
two effects, namely expansion of hot nucleus ( leading to 
a decrease of the absolute values of chemical potentials 
compared to the zero temperature values as indicated by 
a solid curve ) and chemical separation into an isospin 
symmetric heavy fraction and an isospin asymmetric nu- 
cleon gas, leading to stronger dependence of [i n — jjb p of the 
nucleon gas ( dilute phase ) on N/Z of the system 0,0 
thus resulting in higher sensitivity of the observed yield 
ratios. Such a conclusion is supported by the inhomoge- 
neous isospin distribution observed for the current data 
|l9j where the isospin asymmetry of light charged parti- 
cles exhibits strong sensitivity on N/Z of the quasipro- 
jectile while the isospin asymmetry of heavier fragments 
is rather insensitive. The estimate of (3'T for homoge- 
neous system was obtained assuming p 2 / 3 -dependence 
and using the estimate of free volume ( eq. 52 of ref. 
[23j| ) for each excitation energy bin from the successful 
DIT/SMM simulations. Since the isoscaling parameter 
depends directly on the free nucleon chemical potentials, 
the turning-point at 4 MeV/nucleon can be understood 
as a signal of the onset of chemical separation which re- 
verts the decrease of the free nucleon chemical potential 
consistent with expansion of the homogeneous system. 
This is also consistent with the recent theoretical study 
of thermodynamical properties of finite nuclei of Sanzhur 
et al. where the hot nucleus ( at a given initial excita- 
tion energy ) can expand up to the turning point where 
the central part begins to contract again while the in- 
creasing number of nucleons form the dilute component 
surrounding the dense central part. The values of the 
temperature characterizing the turning point in the cen- 
tral density for a given excitation energy reproduce well 
the caloric curve [2(| for the data investigated here. The 
existence of the turning-point explains the difference of 
caloric curves ( around 4 MeV/nucleon ) obtained in p| 
using double isotope ratio and mirror nucleus tempera- 
tures, former free of assumptions on chemical potentials ( 
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the same applying to Fig. ) while latter using the ap- 
proximation [ip. n w —Sp. n ( thus artificially lowering the 
values of (3'T around the turning-point ) . When estimat- 
ing the density of the homogeneous system at the turning 
point ( assuming /9 2 / 3 -dependence of (3'T ), it is consis- 
tent with theoretical estimates of the position and shape 
of the spinodal contour, typically at almost constant to- 
tal density 0.6po in the wide range of asymmetries 0. 
Thus, inside the spinodal region the homogeneous nu- 
clear medium is quickly replaced by the inhomogeneous 
system. 
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In summary, the isoscaling is investigated using the 
fragment yield data from fully reconstructed quasi- 
projectiles observed in peripheral collisions of 28 Si with 
I24,ii2g n a £ p ro j ec tile energies 30 and 50 MeV/nucleon. 
The excitation energy dependence of the isoscaling pa- 
rameter (3' is observed which is independent of the beam 
energy. For a given reconstructed quasi-projectile pro- 
duced in reactions with different targets no isoscaling is 
observed. Thus, the magnitude of the isoscaling is de- 
termined by the difference of initial quasi-projectile sam- 
ples produced in reactions with different targets leading 
to different level of N/Z-equilibration. The excitation en- 
ergy dependence of the observable (3'T, corrected for the 
1 /T temperature dependence, exhibits a turning-point 
at 4 MeV/nucleon which can be interpreted as a signal 
of the onset of separation into isospin asymmetric dilute 
and isospin symmetric dense phase. The initial trend of 
expansion of the homogeneous system is reverted by the 
transition into an inhomogeneous system with two phases 
having different proton and neutron concentrations. A 
novel experimental information on the equation of state 
of the asymmetric nuclear matter is thus obtained. 



FIG. 4: The dependence of observable (3'T on relative exci- 
tation energy of the quasiprojectile. Open squares and circles 
show experimental data obtained using the double isotope ra- 
tio (a) and Y( 3 H)/Y( 3 He) slope (b) temperature at 50 and 30 
MeV/nucleon, respectively. Full squares represent combined 
data from both energies. Full and dashed horizontal lines 
show estimate of (3'T obtained when assuming fi Pt „ ~ —S p . n 
at 50 and 30 MeV/nucleon, respectively. Solid curve repre- 
sents the expected trend for homogeneous system. For details 
see text. 
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